2005 International Linear Collider Workshop - Stanford, U.S.A. 



FERMILAB-Conf-05-318-T 
DESY 05-122 
WUE-ITP-2005-007 



Indirect Determination of the Higgs Mass Through Electroweak Radiative 
Corrections 

A. Freitas 

Fermi National Accelerator Laboratory, Batavia, IL 60510-500, USA 
M. Awramik 

DESY, Platanenallee 6, D-15738 Zeuthen, Germany; 

Institute of Nuclear Physics PAS, Radzikowskiego 152, PL-31 342 Cracow, Poland 
"j^ ; M. Czakon 

' Institut fur Tfieoretische Physik und Astrophysik, Universitat Wurzburg, Am Hubland, D-97074 Wurzburg, Germany; 

. Institute of Physics, University of Silesia, Uniwersytecka 4, PL-40007 Katowice, Poland 



Electroweak precision observables allow stringent tests of the Standard Model at the quantum level and imply inter- 
esting bounds on the mass of the Higgs boson through higher-order loop effects. Very significant constraints come 
especially from the determination of the mass of the W boson and from the effective leptonic weak mixing angle. After 
shortly reviewing the status of theoretical computations of the W mass, the new calculation of two-loop corrections 
^ , with closed formion loops to the effective leptonic weak mixing angle is discussed in detail. The phenomenological 

implications of the new result are analyzed including an estimate of remaining uncertainties. 

O 1. INTRODUCTION 

O 



In recent years, the Standard Model of electroweak interactions has been confirmed experimentally with outstanding 



rS i' success. Not only was almost all the particle content discovered at accelerator experiments, but their properties and 

^ interactions have been measured with high precision, in agreement with the model prediction. The only missing piece 

D ■ is the Higgs boson, which is responsible for electroweak symmetry breaking. However, even today we can obtain 

'■^ [ meaningful constraints on the Higgs boson mass from electroweak precision measurements. Due to the impressive 

. ^ ' accuracy of some of these experimental results, they are sensitive to electroweak radiative corrections at the next- 

■ to-leading (NLO) and sometimes ncxt-to-next-to-leading (NNLO) level, and thus depend on the impact of the Higgs 
boson entering in the loops. 



Two of the most important quantities in this respect are the mass of the W boson, Mw and the sine of the leptonic 
effective weak mixing angle sin^ ^efj** ■ The H^-boson mass can be inferred from the muon decay constant , which 
is generated through virtual W-boson exchange, so that cx The effective weak mixing angle, on the other 

hand, reflects the ratio of the vector and axial- vector couplings, Vf and a/, of the Z boson to fermions (/) at the Z 
boson pole: 



(1) 



Since these couplings can be measured most precisely for leptons, the leptonic effective weak mixing angle sin^ 9^^ 
is usually taken as a reference. 

The current experimental world average for the M^-boson mass is Mw = (80.425 ± 0.034) GeV [1]. Recently, a lot 
of progress has been made towards establishing accurate theoretical prediction for Afw The best result [2] includes 
the complete two- loop corrections [3, 4] and some three- loop contributions [5, 6]. The remaining theoretical error 
is estimated to be 5My^ ~ 4 MeV, which is well below the current experimental uncertainty. Still, the electroweak 
two-loop corrections total to ^ 30 McV and are thus mandatory for electroweak precision analyses. 
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The effective weak mixing angle sin^ 0^^* is mainly derived from various asymmetries measured around the Z 
boson peak at e+e~ colliders after subtraction of QED effects. The current experimental accuracy, sin^ ^^g** = 
0.23150±0. 00016 [1] implies strong indirect constraints on the allowed range for the Higgs boson mass Mh- Therefore 
it is important to develop precise theoretical calculations for this quantity. 

Usually, sin^ (^^^ is computed as a function of the electromagnetic coupling a, the muon constant and the 
masses of the Z boson, Mz, and the top quark, mt (other fermion masses are numerically irrelevant). As explained 
before, A^fw is calculated from G^, but in addition to these corrections, the computation of sin^ f?^*^* involves the 
corrections to the Z vertex form factors. The latter are expressed by the quantity k = 1 + Ak, in such a way that 
the effective weak mixing angle can also be written as: 

sin^ Q'X = (1 - M^/M|) (1 + Ak) , (2) 

and at tree-level, Ak = 0. Higher-order corrections to sin^ have been under extensive theoretical study over the 
last two decades. Besides the one-loop result [7, 8], two- and three-loop QCD corrections are available [5, 9, 10], but 
for the clcctrowcak two-loop contributions only partial results were known. By means of a large mass expansion in 
the heavy top quark mass, the formally leading 0{o?mX) [11, 12] and next-to-leading 0{o?m\M'^ [13] terms were 
computed. A part of the missing two-loop contributions was incorporated by the complete electroweak two-loop 
corrections to Afw [3, 4]. While these corrections effected a shift in Mw of 4 MeV compared to the previously known 
0(o?rn^M'^ contributions, the induced shift in sin^ was very sizable, (5 sin^ 0^^* = 8 x 10~^, thus implying that 
the missing two-loop terms in the form factor Ak can be of similar order. 



2. ELECTROWEAK TWO-LOOP CORRECTIONS TO sin^ B^X 

As a first step towards completing the electroweak two-loop corrections to sin^ ^Ifr**: results for the fermionic (i.e. 

diagrams with closed fermion loops) two- loop corrections were presented recently [14]. The genuine two- loop vertex 
diagrams are represented by the generic topologies in Fig. 1. Higher-order corrections to the process e+e^ ff 
near the Z pole can be consistently computed by performing an expansion of the amplitude around the complex pole 
A^l = M| - iMzTz, 



R 



A[e+e- ^ff] = j-j^ + S + {s~Ml)S' + ... (3) 



Here Fz is the Z decay width. After subtracting contributions from s-channel photon exchange and 'y-Z interference, 
the vertex corrections form factor at NNLO is derived to be 



(2) _ af_ 

'*f . fnis„ . fni (o)x 



(4) 



where the superscripts in parentheses indicate the loop order. In this quantity, IR- divergencies from QED contri- 
butions drop out, which involves a delicate interplay between one- and two-loop terms in the form factors. The 
UV-divergencies are cancelled by on-shell renormalization. The relevant counterterms are derived using the methods 
of Ref. [3]. 
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□ 




Figure 2: (a) General representation of a two-loop scalar diagram with self-energy sub-loop, (b) Reduction of triangle sub-loop 
to self-energy sub-loop by means of Feynman parameters. 



The new part of this work is the computation of the two-loop Z] ] vertex corrections, which are treated with 
two independent technical methods. The first method uses large mass expansions for the diagrams with internal 
top-quark lines and the differential equation method for diagrams with only light fermions f ^ t, the masses of 
which are neglected. Contrary to previous work [13], the expansion in x = M^/m1 is performed to high precision, 
by executing the series to the order a;^°, reaching an overall relative precision of ~ 10~^ of the final result. The 
coefHcjients of the expansion are 2-loop tadpole and 1-loop vertex diagrams, which can be evaluated efficiently with 
well-known analytical formulae. The contributions from diagrams without top-quark propagators involve only two 
independent scales, Mi and Afw, allowing a fully analytical treatment. Even for the limited set of diagrams with 
closed fermions loops, a large number of scalar integrals with non-trivial structures in the numerator are involved. 
They can be reducedj^o a set of scalar master integralsjjy using integration- by-parts identities [15]. Owing to size 
of the linear equation system associated with this reduction procedure, the algorithm has been implemented in the 
dedicated C-|--f library DiaGen/IdSolver [16], which performs the necessary steps in a highly automized way. 
Analytical results for these master integrals are obtained by the differential equation method [17]. This is illustrated 
by the followin^~^xample , 




Here the thick lines represent massive propagators with mass m, the thin lines denote massless propagators and 
p is the momentum flowing into the vertex. D is the dimension of dimensional regularization. The momentum 
derivate of the scalar integrals on the left-hand side results in the same integral and simpler integral topologies on 
the right-hand side. Feeding in analytical results for these simpler integrals, the differential equation can be solved 
in terms of generalized poly logarithms. All integrals were also checked by using low- momentum expansions. 

The second method makes use of numerical integrations based on dispersion relations. A scalar two-loop integral 
with a self-energy sub-loop as in Fig. 2 (a) can be expressed as [18] 

f°° [ \ 1 1 
TN+i{pt;m^) = - dsABo{s,m%,m%^-^) d^q— — 2 ' ' ' T"^ [ ] 2 ' (6) 

Jso J ~ s [q+piY - m{ [q + pi + . . . + pN-i)- m%_^ 

where A.Bq is the discontinuity of the scalar one-loop self-energy function. The second integral can be evaluated into 
a standard A'^-point one-loop function, leaving the integration over s to be performed numerically. In general, one 
can also introduce dispersion relations for triangle sub- loops [19], but it is often technically easier to reduce them to 
self-energy sub- loops by introducing Feynman parameters [20], 

[{q + pi)^ -mj]-^ [(<7 + P2)'-m2]-i dx ^ -^2 _ ,^2]-2 

p — xpi + (1 — x)p2, ffi^ — X nil + (1 ~ a;)TO2 ~ 2;(1 — x){pi — ^2)^- 

This is indicated diagrammatically in Fig. 2 (b). The integration over the Feynman parameters is also performed 
numerically. As a result, all master integrals for the vertex topologies can be evaluated by at most 3-dim. numerical 
integrations. Before performing the numerical integrations, possible UV- and IR-divergencies need to be subtracted 
from the integrals. While this second method is applicable to two-loop vertex corrections with an arbitrary number 



PSN 0610 



«1> 

CN 



0.2325 



0.232 



0.2315 



0.231 




200 



800 



1000 



400 600 
A/h [GeV] 

Figure 3: Prediction for sin^ ^eff** including two-loop corrections compared to the current direct experimental measurement, 
with la bands from experimental input. The chosen input parameters are Mz = (91.1876 ± 0.0021) GeV, Fz ~ 2.4952 GeV, 
mt = 178.0 ± 4.3 GeV, mt = 4.85 GeV, Aa(Af|) = 0.05907 ± 0.00036, a,{Mi) = 0.117 ± 0.002, = 1.16637 x 10"^ GeV'^ 



of mass scales, it is slower and leads to much large expressions than the first method. Nevertheless it provides an 
important check of the result. 

Special caution is needed for the diagrams with a fermion triangle loop (see the third diagram in Fig. 1), which 
involve the 75 matrix. In dimensional regularization, it is not possible to fulfill the two relations {7^^,75} — and 
Tr(7"7''7'''7''75) = Aie"^''^ simultaneously. As in Ref. [3], the contributions resulting in e-tensors were therefore 
evaluated in four dimensions, based on the observation that these terms are free of UV-divergencies. Potential soft 
and collinear divergencies of single diagrams are regulated using a photon mass, with a subsequent careful expansion 
for zero photon masses. 



3. RESULTS 

The new result for the fermionic two- loop corrections is combined with corrections of order 0{a), 0{aas) [9], 
C(aas) [5] and leading three-loop terms of 0{aagmf) and 0{a'^m^) [6]. Reducible terms of the same order are 
taken into account, but no resummations are preformed. The most precise prediction of sin^ ^eff* obtained as 
a function of the muon decay constant G^, from which the M^-boson mass is calculated by including the radiative 
corrections to Mw as given in Ref. [2] . 

Fig. 3 shows the final result for sin^ ^eff * ^ function of the Higgs mass compared to the current experimental 
value. Included in the plot are the error bands due to the uncertainties of the experimental input parameters 
entering into the theoretical prediction and of the direct measurement of sin^^^,*^*. As evident from the figure, the 
identification of computation and measurement for sin^ ^^g** favor relatively small values for A/jj. With the new 
two-loop corrections, the best-fit value for A/h moved from 148 GeV (using formulae of Ref. [13]) to 168 GeV (for 
mt = 178 GeV). 

The numerical result for the leptonic effective weak mixing angle has been published in Ref. [14] as a parametric 
fitting formula that is accurate in the range 10 GeV < A/h < 1 TeV. It has also been implemented into the newest 
version 6.42 of the program Zfitter [21], with some changes recently discussed in Ref. [22], and was used in the 
latest release of electroweak precision global fits of the Standard Model. The impact of the new result for sin^ ^Iff * 
shifts the 95% confidence level upper bound on the Higgs mass upwards by 23 GeV to 260 GeV [1]. 

Together with the inclusion of the new two- loop result in the prediction for sin^ ^eff * j ^ assessment of the uncertain- 
ties from missing higher order contributions is required. Since for practical purposes sin^ ^rff * given as a function 
of the muon decay constant G^, it is useful to evaluate the theoretical error for this parametrization, i.e. combining 
the radiative corrections to AZ-yv and the Z vertex. A simple method to estimate the higher order uncertainties 
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assumes that the perturbation series follows roughly a geometric progression. This presumption implies relations like 
©(a^cts) = 0(0^) /0{a) 0{aas). With this method one obtains the following errors for Mh between 10 and 1000 
GeV in units of 10~^: between 2.3 and 2.0 for the 0{a'^as) contributions beyond the leading term, between 1.8 
and 2.5 for 0{a^), between 1.1 and 1.0 for 0{aa^) and between 1.7 and 2.4 for ©(a^ttg). The missing bosonic 0{a^) 
corrections cannot be appraised from geometric progression. However, considering they have a prefactor but no 
specific enhancement factor, they arc estimated to be about 1.2 x 10~^. To account for possible deviations from the 
geometric series behavior, an overall factor \/2 was included to arrive at a total error of (5thsin^ O^^g*^ ~ 4.9 x 10~^. 

Alternatively, the error from a higher-order QCD loop can be assessed by varying the scale of the strong coupling 
constant as or the top-quark mass mt in the MS scheme in the highest available perturbation order. The scale 
variation leads to an error estimate of 0.1 to 3.9 x 10~^ for the ©(a^as) corrections and of less than 10~^ for the 
0{aal) contributions. These numbers are of the same order as the estimated errors from the geometric progression 
method, so that the total error given above seems to be fairly reliable. 

The new error estimate was used in the latest electroweak global fits [1] and lead to a reduction of the width of 
the well-known blue band, which indicates the theoretical error in the indirect determination of the Higgs mass. 

4. CONCLUSIONS 

In this contribution, recent progress in the calculation of higher-order corrections to the most important electroweak 
precision observables and their impact on the indirect determination of the Higgs mass was reported. 

The complete fermionic O(a^) corrections to the leptonic effective weak mixing angle sin^ 6^^q^ have been calculated 
and numerical results were presented. As an additional check, the computation of the two-loop vertex integrals was 
performed with two independent methods. The new result, together with an estimate of the remaining theoretical 
error, was included in the latest version of the program Zfitter and used for the latest global electroweak fits of 
the Standard Model. 

The calculation of the remaining bosonic electroweak two-loop corrections is currently in progress and will be 
available soon [23]. Furthermore, we are working on adapting the new results for quark final states, where particular 
attention has to be paid to the Zbb vertex, since it includes additional massive top-quark propagators [23] . 
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